New Physics contributions to the lifetime difference in D°-D mixing 
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The first general analysis of New Physics contributions to the D°-D lifetime difference (equiv- 
alently AFd) is presented. The extent to which New Physics (NP) contributions to |AC| = 1 
processes can produce effects in AFd, even if such NP contributions are undetectable in the current 
round of D° decay experiments, is studied. New Physics models which do and do not dominate the 
lifetime difference in the flavor SU(3) limit are identified. Specific examples are provided. 

PACS numbers: 



Quantum mechanical meson-antimeson oscillations are 
sensitive to heavy degrees of freedom which propagate in 
the underlying mixing amplitudes. The observation of 
mixing in the K° and Bd systems thus implied the ex- 
istence respectively of the charm and top quarks before 
these particles were discovered. In like manner, by com- 
paring observed meson mixing with predictions of the 
Standard Model (SM) modern experimental studies have 
been able to constrain models of New Physics (NP). 

Which system of mixed mesons is likely to produce ev- 
idence for NP? It has become clear from B-factories and 
the Tevatron collider that hopes for spectacular NP con- 
tributions to Bd and B s oscillations have come to naught. 
The large SM mixing succesfully describes all available 
experimental data. The only flavor oscillation not yet 
observed is that of the charmed meson D°, where SM 
mixing is very small and the NP component can stand 
out. Q 

Charm mixing arises from |AC| — 2 interactions that 
generate off-diagonal terms in the neutral D mass matrix. 
To second order, the D°-to-D matrix element is 
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where Tt w : is the effective |AC| = 1,2 hamiltonian. 
The most natural place for NP to affect mixing ampli- 
tudes is in the |AC| = 2 piece, which corresponds to a 
local interaction at the charm quark mass scale. This 
local interaction cannot, however, affect ATd because it 
does not have an absorptive part. 

Let us introduce standard notation for ATd and AMd 
by employing the dimensionless forms, 
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Given CP-conservation, we can express y as an absorptive 




FIG. 1: Loop diagram for D° -> D 



part of Eq. (TI| , 
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where p n is a phase space function that corresponds to 
charmless intermediate state n. This^rclation shows that 
APd is driven by transitions D°,D — > n, i.e. physics 
of the |AC| = 1 sector. It turns out that experimentally 
observed D° decays agree reasonably well with SM es- 
timates yL To date, no clear signals of NP have been 
observed |4|. As such, it is currently accepted that ATd 
is dominated by the SM contribution. In this Letter, we 
show that this is not necessarily so and consider several 
NP models to illustrate our point. 

Consider a D° decay amplitude which includes a small 
NP contribution, A[D° -> n] = A ( n M) + A n NP} . Here, 
j4.£ NP ' ) is assumed to be smaller than the current exper- 
imental uncertainties on those decay rates. Then it is a 
good approximation to write Eq. in the form 



y- 



y^ A (SM) A (SU) +2 yPn A (NV) A (SM) (4) 



r D 



The first term in this equation corresponds to SM inter- 
actions at both vertices in Fig. ^ whereas for the second 
term, there is one SM vertex and one NP vertex. 

The SM contribution to y is known to vanish in the 
limit of exact flavor SU(3) [§(• Moreover as was shown 
in 0, the first order correction is also absent, so the SM 
contribution arises only as a second order effect. Thus, 
those NP contributions which do not vanish in the flavor 
SU(3) limit must determine the lifetime difference there, 
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even if their contibutions are tiny in the individual decay 
amplitudes. The same reasoning can be applied to x since 
a dispersion relation relates x to y • 

Of course, flavor SU (3) symmetry is broken in the real 
world. Just how large this effect is on D°-D° mixing in 
the SM is controversial, with estimates for y ranging from 
a percent to orders of magnitude smaller 13 The 
current experimental bounds on y and x are |3J 

y < 0.008 ± 0.005 , x < 0.029 (95% C.L.) . 

A NP | AC | = 1 interaction can have a measureable effect 
on the value of y (and of x) if the true SM value for y 
does not near the top of the range of predictions. We 
shall assume that this is the case. 

Using the completeness relation and Eq. ®, the NP 
contribution to the D°-D lifetime difference becomes 
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We represent the NP AC = — 1 hamiltonian as 1 

Knp^ 1 = E D W PiMQi + C 2 (^)Q 2 ] , (6) 
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Qi = uJig, q J T 2 c i , Q 2 = wJig, qjT 2 Cj , 

where the spin matrices r^a can have arbitrary Dirac 
structure, Ci j2 (m) are Wilson coefficients evaluated at en- 
ergy scale fi and the flavor sums on q, q' extend over the 
d, s quarks. We shall expand the time-ordered product 
of Eq. (|SJ in an operator product expansion (OPE), i.e. 
in terms of local operators of increasing dimension 9] . 

The leading term in the OPE is simply that depicted 
in Fig.^ For a generic NP interaction, we calculate that 
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where {K a } are combinations of Wilson coefficients, 
K x = (CxCxN c + (dC 2 + CxC 2 )) , K 2 = C 2 C 2 (8) 

with the number of colors N c = 3. The operators {Of } 
in Eq. Q) are defined as 
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1 Throughout this paper, we reserve indices i,j,k,£ for color. 



o 



O 



ijki _ — 



O 



ijki 



O 



ijki 



u k T f _ l/ i) c r 2 c j uiTiT 11 ^ 



(9) 



where p c is the charm-quark momentum operator, 
1ij.Pl, Pl = (1 + 75)/2 and later we also use Pr = (1 
75)/2. The coefficients I a (x,x') in Eq. J7J are 
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where k* = (m c /2)[l - 2(x + x') + (x - x') 2 } 1/2 with 
x = m 2 /m 2 and x' = m 2 , /m 2 . 

Equations l(7|l- i|10[l represent the basic formulas in our 
analysis. Hereafter, we take uid = and express re- 
sults in terms of the Wolfenstein parameter A = ~V US = 
-V cd ~ 0.22 and x s = m 2 Jm 2 c ~ 0.006. All our results 
are presented to leading order in x s . Finally, predictions 
using NP masses and couplings other than the ones as- 
sumed here can be obtained via simple scaling. 
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The SM hamiltonian H^m 1 ^ s recovered in Eq. © 
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T M , and the known SM result j/g M ' easily follows, 
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Note that this contribution is suppressd by six powers 
of m s and is therefore tiny, y^^^ ~ 10~ 8 1]. This 
is because the GIM mechanism requires four chirality 
flips (strange quark mass insertions) for the intermedi- 
ate quarks plus additional hclicity flip due to the pseu- 
doscalar initial state. 0. 

Examples of New Physics Models 

In what follows, we distinguish between NP models 
which vanish in the limit of SU(3) flavor symmetry from 
those which do not. 

Nonzero SU(3) Limit: For NP models with flavor- 
dependent couplings D qq i , it is possible to obtain contri- 
butions that are nonzero in the flavor SU(3) limit. For 
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these, the main contibution will come from the opera- 
tors 0i,2 (as 03,4,5 are suppressed by powers of m s /m c ). 
The two most common scenarios involve (V-A)<£>(V-A) 
and (S-P)<g>(S+P) couplings. 

As an example, consider a NP model whose low energy 
effective hamiltonian is represented by a four-fermion op- 
erator with vectorial left-handed interactions, i.e. T\ = 
F 2 = I^Pl and D qq > = X qq >/A 2 , where A is the NP mass 
scale. We find 
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where C = V2G F m 2 c /(3TrM u T D ) and 

A = X s d — X (Xdd — A ss ) — A 2 (Xds + X s d) 
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is the combination of NP couplings to the s, d quarks. 

It follows from Eq. I|13fl that if all NP couplings X qq i 
are of the same order, ?/vlh is nonzero in the flavor SU (3) 
limit. The result for a penguin-like NP contribution can 
be obtained if one sets X s d = X ds = 0. In this case the 
same conclusion holds if Add ^ X ss (which is however not 
easy to arrange) unless a (generally tiny) |>l£ NP ' ) | 2 term 
is also included in Eq. Q. In what follows, we will be 
neglecting QCD running of the local operators generated 
by the NP interaction (i.e. C\ = and C 2 = 1). 

Models with extra vector-like quarks: Consider a model 
of the above type which extends the SM by including new 
singlet quarks in a vector-like representation |l0| . In this 
instance, the Z-boson has additional flavor-changing cou- 
plings. For example, assume both up- type and down- 
type exotic quarks U a ,i, D a i are present (indices a,i 
denote flavor and color respectively). Then the flavor- 
changing couplings are described by 
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with flavor-changing couplings in both up and down sec- 
tors. The lifetime difference for this model can be ob- 
tained from Eq. I|12|) by substituting X s d = ui„U^ , 
Ads = Add = A ss = and A = ^2/Gp- This model 
is well-constrained from measurements of the mass dif- 



ferences in KK and DD mixing. For U 
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nitude as y§ M above. It is worth noting that the little 
Higgs models, which have similar low-energy signatures, 
are not constrained by the measurements of lifetime dif- 
ference, as they do not have flavor-changing couplings for 
the down quark sector (flavor-conserving contributions 
cancel out in y). 

SUSY without R-parity (slepton exchange): Another 
example of a contribution which survives in the flavor 
SU(3) limit is SUSY without R-parity 11]. In this 



model, there are flavor-changing interactions of slcptons 
that can be obtained from the lagrangian 
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as well as the interactions mediated by squarks discussed 
below. The slepton-mediated interaction is not sup- 
pressed in the flavor SU(3) limit and leads to 
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mass, A is given by Eq. i|13|) with X s d = X' il2 X' i21 < 
1 x 10- 9 , A ss = X' m X' m < 5 x HT 5 , A^ = X' l21 X' l22 < 
5 x 10" 



, Ads = A' 2ll ^ 22 < 5 x 10- 2 11], and (Q>) is 
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Operators with a tilde are obtained by swapping color 
indices in the charm quark operators. Using factoriza- 
tion to estimate matrix elements of the above operators 
and taking for definiteness Mj= 100 GeV, we arrive at 
ypt — —3.7%. The contribution due to squark exchange 
vanishes in the flavor SU(3) limit and is given below. 

Zero SU(3) Limit: There are several reasons that 
some NP models vanish in the flavor SU(3) limit. First, 
the structure of the NP interaction might simply mimic 
the one of the SM. Effects like that can occur in some 
models with extra space dimensions. Second, the chiral 
structure of a low-energy effective lagrangian in a par- 
ticular NP model could be such that the leading, mass- 
independent contribution vanishes exactly, as in a left- 
right model (LRM) . Third, the NP coupling might explic- 
itly depend on the quark mass, as in a model with mul- 
tiple Higgs doublets. There, the charged Higgs couplings 
explicitly depend on quark mass. However, most of these 
models feature second order S'?7(3)-breaking already at 
leading order in the l/m c expansion. This should be 
contrasted with the SM, where the leading order is sup- 
pressed by six powers of m s and the second order only 
appears as a l/rn^-order correction. 

Left-right models: Left-right models (LRM) provide 
new tree-level contributions mediated by right-handed 
(W (R) ) bosons [13. The relevant effective lagrangian is 
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where Vk are the coefficients of the right-handed CKM 
matrix. This leads to a local AC = — 1 hamiltonian as 
in Eq. © with Ti = T 2 = j^Pr. Since current experi- 
mental limits allow masses as low as a TeV [3j, a 
sizable contribution to y is quite possible. Using Eq. J7J, 
we obtain 
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where Clr 

g^/8M$y , C 12 are the SM Wilson coefficients and the 
operators appearing in Eq. (|19|) are given in Eq. (|f 7J) . 
Using we obtain numerical values for two possi- 
ble realizations: (i) 'Manifest LR' (V^ L ) = V^) gives 
y LR = _4.8 • lCT 6 with M Wr = 1.6 TeV and (ii) 'Non- 



manifest LR' (V, 
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-8.8 • 10" 5 with 



gives y LR 

Mw R = 0.8 TeV. In both cases we take g-R, — </l- 

Multi-Higgs models: A popular realization of this type 
is the two Higgs doublet model (2HDM) with natural fla- 
vor conservation. This model provides new tree-level con- 
tributions mediated by charged Hig gs b osons and leads 
to the local four fermion interaction |l3| 
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where the vertices Ti^ are 

Ti=m ql cot V uq > P R - m u tan/3 V uq > P L , 

T 2 =m q cot (3 V* q P L - m c tan V* q P R . (21) 

There are four possible contributions involving the var- 
ious terms in V\.2- However, three of these, including 
the potentially large tan 2 /3 term, vanish for assorted rea- 
sons {e.g. flavor cancellation, zero matrix element). This 
leaves 

ychH = ^P^ff cot 2 (3 [d + C 2 ] (Q) , (22) 

where (Q) is as in Eq. Ill Assuming values Mr = 
85 GeV and cot/3 = 0.05, consistent with constraints 
obtained from the observation of B — > tv t |3J] , we obtain 
y chH o± 2 ■ 10- 10 . 

SUSY without R-parity (squark exchange): The 
baryon-number violating squark exchanges arise from the 
lagrangian [ll| 
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This interaction has the same Dirac structure as the LRM 
discussed earlier and leads to 
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where C" = Gf A/ (2mDirT), M~ is a squark mass and 

sq 

the matrix elements (Q 1 ), (Q 1 ) are given earlier. Us- 
ing factorization for the matrix elements, A^A'^ ~ 
3 • 10^ 4 and taking M~ = 100 GeV, we arrive at the 
result v' a ~ 6.4 • 10" 6 . 



In conclusion, we have explored how NP contributions 
can influence the lifetime difference y in the charm sys- 
tem. We argued that the NP signal is dominant in the 



formal flavor SU(3) limit. We also showed that, for some 
NP models, it is possible that small NP contributions to 
| AC | = 1 processes produce substantial effects in the 

D°D lifetime difference, even if such contributions are 
currently undetectable in the experimental analyses of 
charmed meson decays. Coupled with a known difficulty 
in computing SM contributions to D-meson decay am- 
plitudes, it might be advantageous to use experimental 
constraints on y in order to test various NP scenarios due 
to better theoretical control over the NP contribution and 
SU(3) suppression of the SM amplitude. 

We thank X. Tata and J. Hewett for clarifying discus- 
sions. The work of E.G. was supported in part by the 
U.S. National Science Foundation under Grants PHY- 
0244801 and PHY-0555304. S.P. was supported by the 
U.S. Department of Energy under Contract DE-FG02- 
04ER41291. A.P. was supported in part by the U.S. 
National Science Foundation CAREER Award PHY 
0547794, and by the U.S. Department of Energy under 
Contract DE-FG02-96ER41005. 



[2] 

[3] 

[4] 
[5] 

[6] 

[7] 
[8] 

[9] 

[10] 
[11] 



[12] 
[13] 

[14] 



Here we use the phase conventions of E. Golowich and 
A. A. Petrov, Phys. Lett. B 625, 53 (2005). 
See F. Buccella et al, Phys. Rev. D 51, 3478 (1995). 
W. M. Yao et al. [Particle Data Group], J. Phys. G 33 
(2006) 1. 

S. Bergmann and Y. Nir, JHEP 9909, 031 
(1999); Y. Grossman , A. L. Kagan and Y. Nir, 
arXiv:hep-ph/0609178 

J. F. Donoghue et al, Phys. Rev. D 33, 179 (1986); 
L. Wolfenstein, Phys. Lett. B 164, 170 (1985); 
E. Golowich and A. A. Petrov, Phys. Lett. B427, 172 
(1998); A. Datta, D. Kumbhakar, Z. Phys. C27, 515 

(1985) ; A. A. Petrov, Phys. Rev. D56, 1685 (1997). 

A. F. Falk, Y. Grossman, Z. Ligeti and A. A. Petrov, 
Phys. Rev. D 65, 054034 (2002). 

A. F. Falk, Y. Grossman, Z. Ligeti, Y. Nir and 
A. A. Petrov, Phys. Rev. D 69, 114021 (2004). 
A. A. Petrov, arXiv:hep-ph/0311371 G. Burdman and 
I. Shipsey, Ann. Rev. Nucl. Part. Sci. 53, 431 (2003); 
S. Bianco, et. al, Riv. Nuovo Cim. 26N7, 1 (2003). 

H. Georgi, Phys. Lett. B297, 353 (1992); T. Ohl, G. Ric- 
ciardi and E. Simmons, Nucl. Phys. B403, 605 (1993); 

I. Bigi and N. Uraltsev, Nucl. Phys. B 592, 92 (2001). 
E.g., see G. C. Branco, P. A. Parada and M. N. Rebelo, 
Phys. Rev. D 52, 4217 (1995). 

E.g., see H. Baer and X.Tata, Weak Scale Supersymme- 
try: From Super-fields to Scattering Events (Cambridge 
University Press, Cambridge, England 2006). 
R. Mohapatra and J. Pati, Phys. Rev. D 11, 566 (1975). 
E. Golowich and T. Yang, Phys. Lett. B 80, 245 (1979); 
J. F. Donoghue and L. Li, Phys. Rev. D 19, 945 (1979). 
R. Barbieri and A. Masiero, Nucl. Phys. B 267, 679 

(1986) . 



